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Abstract: Motivated by the properties of early universe scenarios that produce obser- 
vationally large local non-Gaussianity, we perform N-body simulations with non-Gaussian 
initial conditions from a generalized local ansatz. The bispectra are schematically of the 
local shape, but with scale-dependent amplitude. We find that in such cases the size of 
the non-Gaussian correction to the bias of small and large mass objects depends on the 
amplitude of non-Gaussianity roughly on the scale of the object. In addition, some forms 
of the generalized bispectrum alter the scale dependence of the non-Gaussian term in the 
bias by a fractional power of k. These features may allow significant observational con- 
straints on the particle physics origin of any observed local non-Gaussianity, distinguishing 
between scenarios where a single field or multiple fields contribute to the curvature fluc- 
tuations. While analytic predictions for the non-Gaussian bias agree qualitatively with 
the simulations, we find numerically a stronger observational signal than expected. This 
suggests that a more precise understanding of halo formation is needed to fully explain the 
consequences of primordial non-Gaussianity. 
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1. Motivation 

Non-Gaussianity that originates from the inflationary epoch leaves distinct signatures in 
present-day astrophysical measurements, and therefore provides a unique link to the early 
universe. Interactions of the field(s) sourcing the primordial curvature fluctuations intro- 
duce non-Gaussian imprints in the statistics of the temperature fluctuations in the Cosmic 
Microwave Background (CMB) and of the density fluctuations that collapse into bound ob- 
jects. These effects give us many independent probes of the signals of inflationary physics 
at different redshifts, sensitive to a range of scales. While current measurements from the 
CMB confirm that the spectrum of primordial fluctuations is Gaussian to a remarkable 
part in 10 3 , that bound is still four orders of magnitude away from testing primordial 
non-Gaussianity at the level predicted by slow-roll inflation and more than one order of 
magnitude above the level expected from non- linear post-inflationary processing of the fluc- 
tuations (see eg [1] for a recent calculation). In addition, any deviation from the simplest 
single field slow-roll inflationary scenario, including multiple fields, derivative interactions, 
features in the potential, or non-Bunch-Davies initial conditions (see [2] for a summary) can 
lead to observable non-Gaussianity at levels within current constraints but well above the 
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slow-roll prediction. Upcoming data from the Planck satellite [3] and a variety of galaxy 
surveys [4, 5, 6, 7, 8, 9] have the potential to achieve accuracy on non-Gaussianity at the 
level expected from non-linear evolution alone. For recent reviews, see [10, 11, 12, 13, 14]. 

Primordial non-Gaussianity is most effectively constrained by complementary measure- 
ments from the CMB and Large Scale Structure (LSS). The galaxy power spectrum and 
bispectrum and cluster number counts provide independent statistics with different system- 
atics, sensitive to different qualitative features of the primordial non-Gausianity. Combined 
with the CMB, these observations will constrain a wide range of qualitative features of any 
observed non-Gaussianity (including amplitude, shape, sign and scale dependence) which 
can rule out large classes of inflationary models. 

Non-Gaussianity of the local type (with bispectrum maximum in the squeezed limit; 
k\ pa &2 3> k%) has recently generated a good deal of interest in part because it will be 
especially well-constrained by LSS observations [15]. Even in the case of Gaussian fluctu- 
ations, the statistics of collapsed objects are different from those of the underlying density 
field, and the ratio of the clustering of the two is known as the halo bias [16]. The partic- 
ular coupling of long and short wavelength modes in local non-Gaussianity introduces an 
additional, distinctive correction (proportional to 1/fc 2 ) in the power spectrum of collapsed 
objects which will allow strong observational constraints on local non-Gaussianity [15, 17]. 
From a theoretical point of view, observably large primordial non-Gaussianity of this type 
requires at least two fields to contribute to the scenario - single field inflation alone can 
only generate a bispectrum of the local shape with an extremely small amplitude (of order 
the spectral index of the primordial power spectrum [18]). The complete phenomenology 
of multi-field models is rich but we will show here that there are qualitative differences that 
are observationally distinguishable in the halo bias. We propose here a generalization of 
the local ansatz that is phenomenologically useful and captures the physics of many possi- 
ble multi- field models. The generalized ansatz allows for different types of scale-dependent 
amplitude /nl along with the standard local shape. 

While signatures of primordial non-Gaussianity in LSS can often be predicted analyt- 
ically, accurate comparisons of observables with theoretical predictions require the inter- 
mediate step of numerical simulations to validate or correct any the analytical relations. 
In this paper we build on previous work of Dalai et al. [15] to numerically investigate 
the effect of scale-dependent, local non-Gaussian initial conditions. Interestingly, we find 
theoretically and numerically that the halo bias is sensitive to two different types of scale 
dependence that can constrain and distinguish between inflationary models. However un- 
like in the constant /nl case, the simplest theoretical prediction for the bias in models 
with scale-dependent non-Gaussianity does not fully agree with our numerical results. In 
this paper we will motivate our new non-Gaussian ansatz, present the analytic predictions 
from that model and the associated simulations. We will discuss a possible explanation 
for the discrepancy, which indicates that this problem constitutes an interesting test for 
our understanding of structure growth, although we postpone a detailed analysis for future 
work. 

The paper is organized as follows. In Section 2 we discuss in more detail the motivation 
from inflationary theory. A self-contained and purely phenomenological discussion starts in 
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Section 3, where we use the peak-background split method to demonstrate the qualitative 
ways scale-dependent non-Gaussianity may be observable. We also present forecasts for 
differentiating the bispectra based on the analytic predictions. In Section 4 we show the 
results of numerical simulations, which demonstrate a stronger signal than the analytic 
prediction, and so are encouraging for the observational prospects. We speculate on a 
possible explanation for the discrepancy between theory and simulation and then conclude 
in Section 5. 

2. A bigger family for the local ansatz 

The original "local ansatz" to add non-Gaussianity to the primordial perturbations is 
[19, 20, 21]: 

<D(x) = <Mx) + /nl [*&(x) - (^(x))] + . . . , (2.1) 

where 3>(x) is (minus) the gravitational potential, <3?g( x ) is a Gaussian random field and the 
degree of non-Gaussianity is parameterized by (typically constant) /nl- Here a positive /nl 
leads to a positive skewness in the density perturbations, and so more very large objects, 
in the same sign convention as WMAP [22] . 

In many scenarios the primary effect of the non-Gaussian correction appears as a non- 
zero bispectrum, defined as 

<$(fci)$(fc 2 )$(fc 3 )> = (2vr) 3 SUhi + k 2 + fca)B*(fci, k 2 , fc 3 ) . (2.2) 

For the local ansatz above the bispectrum is 

B*(fei, k 2 , fc 3 ) = /nl [2 P*(fei)P$(fe 2 ) + 2 perm.] , (2.3) 

where as usual we define 

($(fc!)*(fc2)) = (2vr) 3 dUki + fe2)P*(fci) = (2tt) 3 8 3 D (ki + k 2 ) 27r ^| . (2.4) 

The subscript D distinguishes the Dirac delta function from the density perturbation. 
Assuming the spectral index, n s , has no significant fc-dependence, the dimensionless power 
spectrum is given by A| = Ao(k/ko) ns ~ 1 . 

CMB data (WMAP7) already constrain -10 < /nl < 74 at 95% confidence ([25]; 
see also [26, 27]) and could potentially achieve A/nl ~ few from the Planck satellite 
[22, 28, 29, 30]. The best current constraint from LSS comes from the scale-dependent bias 
induced in the galaxy power spectrum, giving —29 < /nl < 69 at 95% CL ([17]; see also 
[31]). 

While the local ansatz is a useful phenomenological tool, it is only a first step to- 
ward modeling and constraining primordial non-Gaussianity motivated by the fundamental 

We use the convention that /nl is defined in terms of (minus) the gravitational potential early in the 
matter era. We caution that there is another convention (used, for example, in [23, 24]) that defines /nl S 
in terms of the gravitational potential normalized to present day amplitude, which is related to the WMAP 
convention used here by /n£ S = /nl(p(z = oo)/g(z = 0)) (~ 1.36/nl in the WMAP7 cosmology), where 
g(z — 0)/.g(oo) is often referred to as the growth suppression factor. 



- 3 - 



physics of inflation. Specifically, the local ansatz resembles the first term in a series that 
arises from the transfer of isocurvature to curvature fluctuations during or at the end of 
inflation. Such a transfer may be due to additional scalar fields during inflation (multi-field 
[32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46]), or after (the curvaton scenario 
[47, 48, 49, 44, 50, 51, 52, 53, 54]), or inhomogeneous reheating [55, 56]. 

There are at least four possible sources of non-Gaussianity that generate bispectra that 
are largely well-captured by the local ansatz shape in the final curvature perturbations. 
First, a spectator field during inflation is not constrained to have a flat potential, so there 
may be intrinsic non-Gaussianity in that field that is not tightly constrained by the slow-roll 
conditions and which can be transferred to non-Gaussianity in the curvature. Second, in 
multi-field models non-linear evolution of curvature modes outside the horizon will generate 
non-Gaussianity in the observed curvature perturbations even if the field(s) themselves have 
no interactions other than gravitational [44]. Third, the conversion of curvaton isocurvature 
fluctuations to curvature after inflation depends on the energy density in the curvaton field, 
which is at least quadratic in the fluctuations and so introduces non-Gaussianity of the 
local type [47]. Finally, loop corrections may, in special cases, generate a scale-dependent 
non-Gaussianity [57]. 

Phenomenologically, we can write a more general ansatz for the bispectrum of (minus) 
the gravitational potential that is factorizable and symmetric in momentum by introducing 
two functions, £ s (k) and £ m (fc): 



This ansatz captures a wide range of physically motivated and perturbatively controlled 
models, where the functions £ S)m are at most weak functions of scale. The notation refers 
to the physical origin of the two functions in inflationary scenarios: £ s , with s for single 
field, is different from one if one of the fields has non-trivial self interactions or nonlinearly 
sources curvature perturbations; £ m , with m for multi- field, is different from one when 
two or more fields both contribute to the power in curvature fluctuations. We will discuss 
several illustrative examples next. 

2.1 Two field inflation 

First, we consider a two field inflation scenario where running non-Gaussianity can be 
obtained, following [34, 43]. In the SN formalism [58, 59], one uses the dependence on the 
number of e- folds of inflation on the fields present to relate the curvature fluctuations to the 
scalar field fluctuations. Even if a field does not source the inflationary Hubble parameter 
H, the point where inflation ends (and so the number of e-folds, N) can still depend on 
the position of the field. Then, we can express the curvature perturbation resulting from 
fluctuations of two fields 4> and a (up to second order) as 



where for simplicity we have assumed one of the non-Gaussian terms (N ttTIT ) dominates the 
other (N^) and = 0. All quantities are evaluated at horizon crossing for the mode 




(2.5) 




(2.6) 
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k, and Nj, is the derivative of the number of e-folds with respect to the field <j>. To gain 
some intuition about the pattern of multiplications and convolutions in this expression, 
recall that in single field inflation the running of = — 5 ~ contributes the 

term proportional to r] to the spectral index. In addition, the form of the quadratic term 
generates the standard result that the bispectrum in the squeezed limit goes like the spectral 
index, n s — 1, evaluated at horizon crossing of the short wavelength (large k) modes [60] 
(although one must take into account pre-horizon crossing non-Gaussianity generated in 
the statistics of the field Sep to get the complete bispectrum correct). In the single field case, 
the amplitude of the non-Gaussianity in the curvature perturbation is small, but it does 
run if the spectral index runs and the dominant term has scale-dependent /nl evaluated 
at the scale of the short modes 2 . 

In a multi-field scenario, the modes for each individual field have fluctuations of order 
H, so that 

<^(k)^(k')> = (Mk)Mk')} = (2vr) 3 ^(k + k')^|| (2.7) 

= (2irf5 3 D (k + k')P(k) . 

where the asterisk is a reminder that H is evaluated at horizon-crossing for each wavenum- 
ber k. Then the total curvature power spectrum can be written 

(C(k)C(k')} = (2vr) 3 ^(k + k') J P ? (fc) (2.8) 

= (2vr) 3 <5 3 ,(k + kOP(A:)(iV i 2 > + iV, 2 T ) 

= (2vr) 3 <5 3 ) (k + k / )[P CW + P C ( CT )]. 
The tree-level bispectrum, assuming (5<j)5a) = 0, is 

B C (k 1> k a ,k 8 ) = \n^) P ^ ] { k 2 f Pdki)Pdk 2 ) + 5 perm., (2.9) 

B$(ki,k 2 ,k 3 ) = £ s (fc3)£m(fci)£m(fc2) P<s>(ki)P$(k 2 ) + 5 perm., 

where N iCT(T depends on non-trivial self (and gravitational) interaction terms of just the 
field a, so we relabel it £ s , with s for single field. The fraction of power in the a field is 
different from one only if both fields contribute significantly to the power in fluctuations 
so we have labeled this function with an m for multi-field. (Otherwise, the bispectrum 
would reduce to the usual single-field expression, where /nl must be of order slow-roll 
- that is, the term quadratic in 6(f) would be most important, giving a bispectrum with 
the same form as the first line of Eq.(2.9) but with the coefficient of the power spectrum 
terms N^^k).) Assuming that the potential $ is defined in the matter era, the precise 
relationship between the first and second lines above is 

~N t<Tff (k) = (2.10) 

p c (k) " uw ' 



2 A detailed discussion of how this is consistent with real-space formulations of the local ansatz can be 
found in [61]. 
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We note that quite generally all the coefficients N^, N t(T , iV CTCr , etc will be scale- 
dependent as the potentials for the fields are not exactly flat. In that sense, in any two-field 
scenario with large local non-Gaussianity, running of the amplitude through the function 
£ s (k) is as natural as running of the spectral index. It may be somewhat fine-tuned to have 
two fields contribute to the amplitude of fluctuations (although this is hard to say in the 
absence of compelling particle physics realizations of inflation), but if they do it is likely 
natural for their potentials to be slightly different so that £ m (&2) is scale-dependent. We 
will parametrize this scale dependence by writing 

(s).(m) 
/ \ n f 

\ k pJ 

where k p is a (theoretically irrelevant) pivot point. 
2.2 Mixed curvaton/inflaton scenario 

Now suppose the curvature perturbation comes partly from a Gaussian inflaton field ((f)) 
and partly from a 'curvaton' field (a) which was a spectator during inflation but contributes 
to the curvature perturbation afterwards [62, 63, 64, 65, 66, 47]. The curvaton naturally 
has a contribution that is quadratic in real space since it contributes proportionally to the 
energy density in its fluctuations. Assuming a purely quadratic potential for the curvaton 
gives 

p a = -m 2 a (a + 8a) 2 5p a = l -m 2 (2a5a + 5a 2 ). (2.12) 

The field fluctuations are still generated during inflation, with amplitude yj (5a 2 ) = H/2-k. 
The quadratic term means that the curvaton can contribute a local-type non-Gaussianity 
with /nl constant and determined by the proportion of energy in the curvaton at the time 
the field decays. 

Then we can write the total curvature field as a sum of contributions from the inflaton 
and curvaton: 

((x) = C*(*) + Ca(x) + pMUx) 2 - (Ux) 2 )) (2.13) 

where the factor of 3/5 enters since f^L is conventionally defined for the matter era po- 
tential (Eq.(2.1)) rather than the primordial curvature. 

Assuming the fields don't couple, the bispectrum takes the familiar local form, but 
now in terms of P(( a )'- 

B c (k u k 2 M) = lfMPc(a)(h)P C (a)(k2) + 5sym] . (2.14) 
If we define the ratio of power contributed by the curvaton 

e(k) - P ^ {k) - M (2 15) 

t [k) - p cw (fc) + p a<p) (k) - p c (k) [2Ab) 
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we can write 

B c (k u k 2 ,k 3 ) = p^(k 1 )^k 2 )P c (k 1 )P < (k 2 ) + 5 perm] (2.16) 
B$(ki,k 2 ,k 3 ) = im{ki)^m{k2)P^{ki)Pq,{k 2 ) + 5 perm 

where we have absorbed f^ L into £ m , that is £ m (k) = \/ /jyx£ (^) ■ Again, we parametrize 

(m) 

the function £ m (k) as a simple power law, £ m (h) oc k n f . 

2.3 Curvaton alone 

If non-Gaussianity comes from the curvaton alone, and a potential other than quadratic is 
considered, the bispectrum can again take the form [67, 68] 

B*(ki,k 2 ,ks) = Uh)P*(h)P$(k 2 ) + 5perm., (2.17) 

where £ s (k) can be parametrized as a power law, at least for some potentials, and 
apparently can have either sign. Inhomogeneous reheating may similarly generate this 
bispectrum [61]. 

2.4 Relation to the spectral index 

The running of £ m (k) is an important physical feature of either type of two-field model: it 
is the evolution of the relative power in the two fields during inflation. Just as the spectral 
index measures the variation of the overall amplitude of fluctuations during inflation, for 
two-field models the bispectral index rvY^ can provide complementary information about 
how the contribution from each field evolves. For some curvaton scenarios there would be a 
link between running non-Gaussianity and large scale power asymmetry in the CMB [69]. 
There is a precise relationship between the spectral index and the bispectral index 



M. 



dlnP c 
din A; 

d^Pc(g) 

dink 



n s -l (2.18) 



dhifm (m) 
TInT""/ =n °~ ns 

Although the running of the bispectrum may have either sign, models with a red tilt for 
the field a are anecdotally more common and in that case we have 

nf ] < -(n s - 1) . (2.19) 

Notice that some of the literature (e.g. [44]) defines /nl(^) = £,m(k) 2 and so quotes 
rif < —2{n s — 1). Here however, we will see that there are two different shifts in the non- 
Gaussian bias, each dependent on one factor of Cm(k), so we define as the running in 
that function. Finally, notice that the spectral index of the observed curvature perturbation 
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depends on the running of both fields <j> and a. If the running of the fields is large enough, 
it will change which field dominates the curvature statistics. 

Whatever the origin of the running in either function £, m ,s, it parametrizes the devia- 
tion from exactly quadratic potentials in either field and so is expected to be generically 
on the order of slow-roll parameters (and should be to avoid substantial corrections to 
this parameterization). We will use somewhat large values of the running to confirm the 
behavior of this type of model in our simulations, but the observational goal should be to 
measure |n^''- m ' > | ~ 0(n s — 1). We discuss the potential of future surveys to reach this goal 
in Section 5. For the standard quadratic curvaton case, nj l > seems more natural (that 

is, non-Gaussianity increases on small scales) while Byrnes et al. [43] found > rij > —0.1 
in a survey of multi-field hybrid inflation models. The sign can be understood if the non- 
Gaussianity is entirely due to non-linear evolution outside the horizon. Then large scale 
modes (which exit earlier) will to be more non-Gaussian than smal scale modes. 3 4 

2.5 A comment on naturalness and completeness 

Given that it is already difficult to convincingly explain one field with a very flat potential, 
we may reasonably ask if the scenarios we are considering are even less likely than the usual 
single-field inflation. It is very hard to answer that question without better fundamental 
models - it may be that where there is one inflaton-like field, there are naturally several 
(especially in higher dimensional models), or not. In inflation, there is a very compelling 
reason why the spectral index should be slightly different from one: old inflation models 
with exact de Sitter space are difficult to connect to the early, hot universe after inflation, 
while slow-roll with the Hubble parameter not exactly constant can have a natural end to 
inflation and a period of reheating. If the slow- roll scenario is right and if two fields are 
present and relevant during inflation, it may be reasonable to expect that they both have 
nearly flat and yet not identical potentials. If one accepts that local type non-Gaussianity 
is natural (or more compellingly, if it is observed), scale dependence is also natural. In the 
absence of a range of compelling high energy models, it is hard to quantify the likelihood 
of any of these scenarios. 

However, from a phenomenological point of view, considering a generalized local ansatz 
is helpful in two ways: first, it argues for a careful analysis of different mass tracers in 
any test for primordial local type non-Gaussianity and second, it provides a test of our 
understanding of structure formation. As we will see, the existing expressions for halo bias 
do not give particularly satisfactory agreement with our simulations. 

The generalized local ansatz above is useful to uncover new observational signatures, 
and it would be interesting to investigate to what extent it holds in more complicated 

3 Previous authors have employed different notation for scale-dependent local models. In particular, 
Byrnes et al, in an extensive discussion of possible multi-field bispectra [61] propose a definition of /nl 
and its running that for our ansatz correspond to /^ raes (fci, fe, ks) = [£ s (A;3)f m (fci)£ m (A;2).P$(fci)P$(A:2) + 
sym]/[P 4 ,(fc 1 )P $ (fc 2 ) +sym] and n^ aes s din |/^f nes (fti = k 2 = k 3 = k)\ / din k = n ( f 3) + 2n ( f m) . 

4 In a discussion of the ability of observations to constrain two-field models of the mixed curvaton/inflaton 
type, Tseliakhovich et al [70] recently defined a variable xi where x\ — and where only the scale- 
independent case was considered. In addition, their function £ is defined differently: £ hcrc = i/(i-|-(£ there ) 2 ). 
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models with more (and coupled) fields. However, even with this ansatz we are still far 
from considering all possible known effects. In most two field models, we expect higher 
order terms (like a £ 3 contribution) to be present in the expression for the non-Gaussian 
curvature. Those corrections are also important for comparison of observation with realistic 
models and have been considered in [71, 72, 73, 74, 75, 51]. In addition, there are other 
possibilities that require something even more general than the symmetric, factorizable 
form. For example, non-gaussianity generated by loop effects can sometimes be large and 
goes like [57] 

(C 3 ) oc f NL (mm{kx, k 2 , k 3 })[P(h)P{k 2 ) + 5 perm.] (2.20) 

In addition, the power-law behavior of our ansatz is a poor model for scenarios with a 
feature at some particular scale, such as [76] (features in the potential) or [77] where non- 
Gaussianity effectively switches on at some scale where a spectator field becomes light. 

Finally, we note that scale-dependent non-Gaussianity may also arise in other ways 
and for other bispectra, but most other examples in standard inflationary models are less 
divergent in the squeezed limit than the local shape is and so have a weaker signal in the 
bias. However, there is a small region of parameter space in ekpyrotic models that seems to 
generate bispectra with scale-dependent amplitudes consistent with current observations, 
and more divergent than the local ansatz [78]. 

3. Generalized local ansatz and large scale structure statistics 

Our ansatz for the factorizable, symmetric, scale-dependent local bispectrum is 

B<t,(kx, k 2 , fe 3 ) = &(fc3)£m(fci)Cm(fc2)-P*(fci)P*(fc2) + 5 perm . (3.1) 
where we parametrize the /c-dependence of the amplitude as 

(s),(m) 

/ k \ f 

£,s,m{k) = £,s,m(kp) yjT~ J ' (3-2) 

(s) (m) 

with \rij | < 1. Ideally, the pivot scale k p can be chosen at a point where the amplitude 
and running of the shape are as close to uncorrelated as possible. We adopt k p = 0.04Mpc _1 
based on analysis for the CMB in [79], although they used a slightly different ansatz for 
the scale dependence. Note, however, that the constraints on £ s ,m(k) will be entirely 
independent of the chosen value of k p . 

Although generically we might expect both functions £ s (fc) and £ m (k) to be present, we 
can consider the two functions separately for simplicity. For that reason, we will compare 
the following two bispectra in what follows: 

B|(fci, k 2 , k 3 ) = Uh)P*(k2)P*{h) + 5 perm (3.3) 

B^(k 1 ,k 2 ,k 3 ) = U{h)U{k2)P^>(k 1 )P^{k 2 ) + 5 perm 

The first line applies to a model where only one field contributes to the curvature per- 
turbations (and the inflationary background is sourced by something else). For example, 
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it is generated by a curvaton model where the potential has terms other than the mass 
term, Eq. (2.17), or from a simplified version of the 5N case from Section 2 (where we 
take £(&), the ratio of power in the two fields, to be constant). Since the curvature pertur- 
bations come only from one field, we label the new function £ s (k) with s for single field. 
Scale-dependence in this function indicates the presence of non-trivial self-interactions (eg, 
deviation of the curvaton potential from exactly quadratic). The second line corresponds 
to a scenario where (at least) two fields contribute to the curvature perturbations, but 
the relevant self-interactions are purely quadratic. For example, this is the mixed infla- 
ton/curvaton model of Eq. (2.16), where the curvaton has only a quadratic potential. The 
label m on £ m (k) indicates that multiple fields contribute to the curvature perturbations. 
Scale-dependence in £ m shows how much the potentials for the fields differ. 

Notice that the first model in Eq.(3.3) has a form that is equivalent to the bispectrum 
one would get from generalizing the local ansatz by 

<&(x) = <D G (x) + /nl * [^(x) - (^(x))] . (3.4) 

This also justifies the single-field label. Scale dependence of this type was studied recently 
in Ref. [80]. 

3.1 Scale-independent non-Gaussianity and bias 

In this section we will work out a prediction for the possible signatures of our generalized 
local ansatz in the halo power spectrum. We will only be concerned with the behavior 
of the power spectrum at very small k, where the deviation from the Gaussian case is 
largest. The matter perturbations 5 at redshift z are related to the perturbations in the 
early matter era potential $ by 

S(ie,z) = M(k,z)$(k) (3.5) 

M{k,z) = ~^- 2 D(z)^-T{k)k\ 
3il m Hfi g{co) 



o 



so that Ps(k,z) = M (k, z)P^(k). Here Qm is the matter density relative to critical, H[ 
is the Hubble constant, D(z) is the linear growth function at redshift z normalized to one 
today, and the growth suppression factor is — 0.76 in the best-fit ACDM model. 

We use the Eisenstein & Hu [81] fit to the transfer function T(k). The variance of density 
fluctuations at redshift z smoothed on a scale R associated to mass M is <J 2 (M, z), defined 
by 



00 dfc. 



a 2 (M,z)=J^ —W R (k) 2 M(k,z) 2 Al(k). (3.6) 

where the power spectrum of the primordial curvature perturbations is given by Eq. (2.4) 
and Wn(k) is the Fourier transform of the top-hat window function. The spatial smoothing 
scale R is related to the smoothing mass scale M via 

M = ^7rR 3 pm,o, (3.7) 

where p m $ is the matter energy density today. We write the combination M(k, z)Wn{k) = 
M R (k,z). 
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3.2 Peak-background split and halo bias 

Halos in N-body simulations are associated with peaks of the initial, linear density field S oc 
k 2 &, whose heights exceed some threshold [82, 83]. The basic idea of the peak-background 
split [84] is to compute the effect of long-wavelength background modes on the heights of 
small-scale peaks, and thereby estimate the large-scale clustering of halos. The procedure 
used in the peak-background split is to perturb a single background mode A3>(fcj) and 
propagate the effect of this perturbation to the height of a peak near threshold. In Gaussian 
cosmologies, where there is no mode coupling, the heights of peaks are simply boosted 
by the density associated with the background mode, AS(ki) oc kfAQ(ki). With non- 
Gaussianity, however, there is mode coupling, so we have to compute how this background 
mode affects shorter-wavelength modes $>(k s ) as well. This clearly involves looking at the 
bispectrum in the squeezed limit -B$(fc;, k s ,—k s — ki ~ — k s ) where ki <C k s . 

Using the argument above, we can predict the consequences of modifying the local 
ansatz to include some form of scale dependence; our discussion here is similar to that in 
[85, 86, 87]. To get a feel for the effect of the scale-dependent non-Gaussianity on the bias, 
notice that we can rewrite the expression for the AN type non-Gaussian field (the first 
line of Eq. (3.3)) in Fourier space as a sum of Gaussian modes <3?g(^) and a non-Gaussian 
piece $_b(A;) designed to recover the single- field model bispectrum: 



where we have dropped the <5f)(k) term which is not important for this discussion. Now 
we can use Eq.(3.8) to consider the effect of some long-wavelength perturbation A<J>(fcj). 
Considering k ~ — q2 ~ k s and qi = ki in that expression we see that 



where in the latter equation, we have used the Poisson equation. We then sum over all 
the short wavelengths to get the total boost in peak height, which then translates into the 
halo excess and halo bias. Since short wavelength modes are essentially modes with length 
scales up to the scale of the object (there is a window function in the integration), we see 
that the presence of the scale-dependent function £ s (fc) implies that collapsed objects have 
a shift in bias with amplitude given by an effective /nl roughly on the scale of the object. 
That is, if £ s (fc) (which in this simple case is like /nl(^) as in Eq.(3.4)) increases on small 
scales, smaller mass objects will have a larger non-Gaussian correction than very massive 
objects. 

For the curvaton type model, the second line of Eq. (3.3), the effect above comes along 
with an additional ^-dependence in the bias: 



$(k) = $ G (k) + $ B (k) 



(3.8) 




A$ B {k s ) = 2&(fc a )A*(fc,)$(fe-) 
AS(k s ) = 2£ s (fcs)A$(fc i )<5(A; s ) 



(3.9) 
(3.10) 



A5{k s ) = 2U{h)U(k s )A$(k l )5(k s ). 



(3.11) 
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In this case, different mass tracers have a non-Gaussian shift in bias, Ab, with an amplitude 
proportional to an "effective /nl" set by their mass (a consequence of £ m (k s )) and a scale 
dependence that goes as k~( 2 ~ n f^ (a consequence of £ m (lq)). 

3.3 Alternative derivation of scale-dependent effects 

The intuitive procedure above leads to essentially the same result as the procedure outlined 
by Grinstein and Wise [88] and further developed by Matarrese et al. [89]. For a generic 
non-Gaussian distribution, Refs. [88, 89] found a way to express the two-point function of 
peaks in terms of a series expansion in correlation functions. Motivated by the results of 
[15], this expression was recently applied to the case of local non-Gaussianity by Matarrese 
and Verde [90] . We can use the same starting point to consider the effects of our generalized 
local ansatz, and express the two-point function for halos of mass M (Ch,M(\xi — x 2 \)) hi 
terms of the n-point functions of the density field smoothed on the associated scale R 



&i,m(|£i-Z2|) 



(2) 



[4 3) (Xl ,X!,X 2 ) Hr> {X1,X2,X 2 )} + ... (3.12) 



where the dots represent higher order terms (both higher correlation functions and higher 
powers of the two- and three-point function) . The collapse threshold 5 C is contained in v = 
5 c /a(M). Now we can use the Fourier transform of the halo auto-correlation, Eq. (3.12), 
to compute the bias: 



PM,h(k,z) 



a(M) 2 (z) 



v*{z) 



a(M) 3 (z) 



P*(k)M R (k,z) 



(3.13) 
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where k = k 2 + k 2 + 1kk\[x and the redshift independent integral is 

1 /"OO 

FR^nfMp) = 87r2(7 ( M )2 / dk ± kfP^MRih) 



x dfi M R (k) 
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P*(fc) l^T 
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k \ f i 



(3.14) 
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+ 



2n 2 a(M) 2 



(m) 

P$(fci) \ k 




(s) (m)\ . 
n, is a 



The second expression is in the small k limit, so k ~ fci. In that limit, J~n(k, n j , 

constant that depends on the smoothing scale. When = 0, this expression reduces 

to that of [90] and is identically one in the small k limit. 

The Lagrangian halo bias b^ for halos of mass M is defined by 

2 



P h (k) = b 2 L Ps(k) = b 2 Lfi 



(3.15) 



where in the second equality the /nl = contribution bi,,o has been explicitly factored out. 
Then from Eq.(3.13) the change in the bias relative to the Gaussian value is 



Ab 



a(M) 2 (z) 



25 r . 



M R (k,z 



~£,s(kp)[£,m(kp)Y 



(m) 



F R {k,n 



(«) „( m )^ 



(3.16) 



In using the definition of Ab above to derive Eq. (3.16) from Eq. (3.12), we have 
expanded the square-root which is not always strictly valid. However, the resulting expres- 
sion agrees with the peak background split (and we find that keeping the square-root yields 
worse agreement between theory and simulation). The results for the constant, small k part 
of the integral in Eq. (3.14), J-n{k) — > J~(M, feci) for the representative two-parameter 
cases are plotted as a function of mass (related to smoothing scale R following Eq. (3.7)) 
in Fig. 1. The functions £ Sjm are normalized to (, s ,m(k p ) = 1 so that the left panel shows 
an effective f^L generated by the scale-dependence for each scenario. The right panel 
compares the prediction for the non-Gaussian correction to the (Lagrangian) bias for the 
single- field and multi- field scenarios. (This label indicates how many fields contribute to 
the curvature perturbations - the inflaton itself may be separate). 

3.4 Summary of analytic results 

We have arrived at the same prediction in both of the previous subsections: a general 
factorizable and symmetric extension of the local ansatz leads to two possible modifications 
of the non-Gaussian halo bias. First, different mass objects may see a different non- 
Gaussian correction that goes roughly like the amplitude of the non-Gaussianity on the 
scale of the object. Second, the power of k appearing in the scale-dependent correction 
can be shifted away from the standard k~ 2 result when there are two fields contributing 
to the curvature power and their relative importance is a function of scale. Either the 
first effect alone or a combination of both may be found, depending on the origin of the 
scale dependence. The most general case has two parameters to characterize the running, 
and one to characterize the amplitude. From the point of view of measurements of bias, 
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M / (h" 1 M s J k (h Mpc" 1 ) 

Figure 1: Left panel: The effective amplitude of the non-Gaussian bias on small scales (/nl) as 
a function of the object's mass for two modifications of the local ansatz. The blue short dashed 
lines are the single field model (only £ s (fc) different from one) and the red long dashed lines are 
the multi-field (only £ m (k) different from one). The upper lines show the effect of non-Gaussianity 
that increases on small scales, with rif = 0.6 or = 0.3 while the lower lines have n f = —0.6 
or n™ = —0.3. All curves are normalized to £ Sim (fc p ) = 1. Right panel: A comparison of the 
correction to the bias of objects of mass 4.4 x 10 14 h~ x M@. The solid black curve is the usual local 

(s) 

ansatz, the blue long dashed curve is the single-field model with = 0.6, the red short dashed 
curve is the multi-field scenario with — 0.3, and the purple dot-dashed curve is the multi-field 
scenario with n . = —0.3. Again, £ s . m (fc p ) = 1. 



these combine into the mass-dependent coefficient of the scale-dependent term, S, and 
the power of k that appears in the denominator. In other words, phenomenologically we 
have (in the small k limit) 

AMfc,M)«^M. (3.17) 

k 2 - n f 

More precisely (and in terms of the Gaussian Eulerian bias b^) 

( m ) 

Ab NG (k,M)-f NL (M,n f ,n f , k p ) ^ J [ c^T(k)D(z)g(0) 



(3.18) 



where 

f^{M,n { f s \nf\ kp) = 6(fc P )K m (^)] 2 ^(^ « l,nf ,n { f) ■ (3.19) 

There is some suggestion, both from simulations and from analytic considerations, 
that there is an additional factor multiplying the expression above for A&ng even in the 
case of constant local non-Gaussianity. For example, Giannantonio and Porciani [91] have 
suggested a multiplication by a factor 
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where Abi is a second order non-Gaussian correction that can be calculated from some 
choice of non-Gaussian mass function (and the subscript / indicates that it is scale- 
independent) 5 . Although it is reasonably well motivated, we do not find that such a 
correction alone substantially improves the fit to our simulations (especially for negative 
/nl)> so we remain agnostic about the analytic form of any additional corrections and 
instead focus on the effects unique to the generalized local form, especially /^f (M). From 
a practical perspective, the coefficient above can be fit from simulation and will not affect 
our conclusions. 



3.5 Forecasts based on the analytic prediction 

We now estimate the ability of future observations to detect slow-roll values of the running 
parameter. Here we present Fisher matrix forecasts based on the analytic predictions 
above. This analysis complements earlier forecasts made for both the scale-independent 
[15, 24, 92, 93], and also scale-dependent [79] models of non-Gaussianity. 

We first consider the two simpler scenarios, given in Eq.(3.3), that isolate the single 
field or multi-field effects and each have only two parameters. These can be obtained by 
setting either of the functions £ s or £ m to one in the general expressions above (and so 
either or ni m ^ is set to zero). The explicit expressions are 



single field : Ab NG (k, M) = ff L (M, nf , k 

ecS m,T As) u n, _ £,s{k p ^ 



2n 2 a{M) 2 J 



3(6g - l)5 c n m H$g(oo) 
c 2 k 2 T(k)D(z)g{0) 
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dk x k'fP^k^Mlih) 



(3.21) 



multi — field 



Ab NG (k,M) = f^ L (M,n i p\k p! 



f cS (M r) (m) h \ - ^ m ^' 
./NlA M > n / > K p) — o_9_/ 



(m) 



3(frg - l)5 c n m H 2 g(oo) 
c 2 k 2 T{k)D(z)g{0) 



2ir 2 a(M) 2 



k, 



We have used only the small k portion of the integral expression from Eq.(3.14) since this 
corresponds to the prediction from the peak-background split, and since the integrand at 
high wavenumbers (e.g. k ~ O(0. 1) /iMpc -1 ) depends on the explicit form of the window 
function. We report constraints on the momentum dependent functions that contribute to 
the integral in fff L (M): 



most general 
single field only 
multi — field only 



fNL{k) = £s(k p )[£, m {kp)] 

k 



(s) . (m) 



(3.22) 



fNL(k) = ^ s (k p ) 



/ 



f N L{k) = [U(k p )]' 



(m) 

k\ n f 



3 We thank Tommaso Giannantonio for detailed correspondence on this point. 
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where the last two lines specialize to the simpler cases of considering only the single-field or 
multi-field effects. Conceptually, f^(M) and its Fourier-space analogue, /nl(&) capture 
the high frequency scale-dependence of non-Gaussianity, k ~ M -1 / 3 . In addition to this, 
the bias has low- frequency scale-dependence oc k f for k <C M ' . This is why 
our expressions for /nl(^) and f^(M) contain only one power of • , even though the 

(m) 

bispectrum has two k n f terms. Fiducial values adopted were /nl(^p) = Cs(k p )^(k p ) = 30 
and rif = while k p = 0.04 Mpc~ as before. We will see in a moment that the true 
best-measured scale from the large-scale clustering of galaxies and clusters is somewhat 
smaller. 

Suppose that we have measurements of the power spectrum using objects (galaxies and 
clusters of galaxies) that have been separated in several mass bins. We assume that the 
covariance matrix of measured Fourier-space overdensities in a given redshift bin centered 
at z is given by 

C ab (k, z) = b(k, M a , z) b(k, M b , z) P(k, z) + S ab (3.23) 

n a (z) 

where the labels a and b refer to mass bins. This equation encodes how to combine ob- 
servations from different mass bins, and also straightforwardly specifies the dependence on 
the parameters of interest /§l(^p) anc ^ n f via Eqs. (3.16) and (3.14). 

The Fisher matrix can now be evaluated in the FKP approximation [94], where infor- 
mation is summed over the redshift bins and wavenumber shells. We have 

f » = a — I Us) * L Tr [c c - c Cj] w (3 - 24) 

where O survey and z max are the solid angle and maximum redshift in the survey respectively, 
V is volume, commas denote derivatives with respect to the non-Gaussian parameters, and 
we have suppressed the dependencies of C on wavenumber and redshift. In practice we 
replace integrals with sums to evaluate this expressions. We neglect the effect of redshift 
uncertainties, but assume thick redshift bins with Az = 0.2. 

For definiteness, we assume a dataset of the quality expected from the Dark Energy 
Survey (DES; [6]), with 

^max — 1 and covering 5000 square degrees; the total volume in this 
survey is about 6.5/i _1 Mpc 3 . We assume k m - m = 0.0001 /iMpc -1 and /c max = 0.1 /iMpc -1 ; 
the latter ensures that all information safely comes from the linear regime. Finally, we 
choose the number density of sources above some mass to correspond exactly to the expecta- 
tion from the Jenkins mass function [95]. Therefore, n a (z) = J^ a ' hiBh (dn/dlnM)(z) dlnM, 
where M 0) i ow and M ai hi g h arid the boundaries of the a-th mass bin. The total number den- 
sity of sources at z = and above 1O 13 ' 5 /i -1 M0 is n ~ 10~ 4 (/iMpc -1 ) -3 . We assume 
a large number of mass bins (forty) in M/ M©, uniformly distributed in log 10 M from 
10 13 - 5 h~ l M & to 10 15 ' 5 /j^M©. 

First considering the single- field case, the error in /nl(&) at any k is given by a simple 
propagation of errors 

1/2 



— — — — Uov ff H j-— L,ov nn + 2— — t— t^— Uov nf 

ofaUkp) dny of NL {kp) q u w 



(3.25) 
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where Cov = F~ l is the covariance matrix of the two non-Gaussian parameters that we 
consider. The errors in /nl(^) are shown in the left panel of Figure 2. 

We are also interested in finding the best-constrained wavenumber, A; unC orr- When 
kp = fcuncorrj then the errors on the parameters £(ftp) and rif'^ m ' are uncorrelated . While 
this best-constrained wavenumber can obviously be read off from Fig. 2, it can also be 
calculated analytically as 

^uncorr = k p exp ( - — — J , (3.27) 

V /NL(fc P )Cov ff y 

where k p = 0.04 /iMpc -1 is the arbitrary pivot in Eq. (3.2). The way that A; uncorr 'runs' 
with changing mass illustrates the point we made in Sec. 3.2 that different mass halos 
probe scale-dependent NG on scales corresponding to those masses. 

We find that the best-constrained wavenumber of our survey, for the single-field model 
and assuming DES-quality data, is ft U ncorr — 0.1 /iMpc" 1 , and the corresponding parameter 
errors at k uncoTr are 

<K/NL(/cuncorr)) ^ 8, a{nf) ~ 0.5 (DES forecast, single - field). (3.28) 

We also find that the error in f-Nh{k uncorr ) is largely insensitive to the fiducial value of 
/nl(^p)) while the error in the spectral index becomes larger for a smaller fiducial 
/nl(^p) (which is expected, since a larger fiducial non-Gaussianity increases the absolute 
change in /nl(& 7^ k uncorr ) for a fixed change in ). 

We repeated the same exercise assuming data of the quality expected from the Large 
Synoptic Survey Telescope (LSST; [9]), with z max = 3 and covering 20,000 square degrees; 
the constraints became 

o-(/NL(fc U ncorr)) ^ 1.7, c^n^) - °- 17 (LSST forecast, single - field). (3.29) 



Next we consider the multi-field model; see the right panel of Fig. 2. As expected, the 
numerical constraints on the amplitude are comparable to the single-field case, however 
the constraints on the running improve, and the best-determined scale moves to a slightly 
lower k: 

^(/NL(^uncorr)) ^ 8, cx(nj. m) ) ~ 0.2 (DES forecast, multi - field) (3.30) 
<K/NL(£:uncorr)) - L7 > <K>4 m) ) - °- 04 ( LSST forecast, multi - field). (3.31) 
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Figure 2: Forecasted constraints on /nl(&) in the model where a single field generates the curva- 
ture perturbations (left panel) and in the multi-field model (right panel). We show forecasts for data 
expected from DES (solid black curve) and LSST (dashed black curve) observations. The wavenum- 
ber at which the constraints are the best is & uncorr , and at this wavenumber the normalization and 
slope of the power law are precisely uncorrelated. The six colored contours on top of each panel 
show the individual constraints from six narrow mass bins uniformly distributed in log 10 M from 
10 13 5 /i _1 Af0 to 10 15 /i _1 M Q (assuming the DES survey). In the single-field scenario, individual 
masses do not break degeneracy between amplitude and running of /nl(&) and only constrain this 
function at a single k value; combined masses are required to break the degeneracy In the multi- 
field scenario, the degeneracy is broken even with halos of a fixed mass. [Note that, in all cases, the 
overall constraints on /nl(^) between different wavenumbers k are strongly correlated, given that 
we are assuming a power law in k.] 

The colored contours in Fig. 2 show the individual constraints from each of the six narrow 
mass bins uniformly distributed in log 10 M from 10 13,5 /i~ 1 Mq to 10 15 }i~ 1 Mq (this is for 
the DES survey scenario and the single- field innaton model) 6 . The thick black curve in 
either panel shows the combined constraint. Note that the combined constraint contains 
the information from the individual bins and the correlations between them (corresponding 
to a 7^ b in Eq. (3.23)). A particularly interesting feature of testing these models with 
primordial non-Gaussianity is that halos of different mass complement in producing the 
overall constraint. For example, inspection of Eqs. (3.16) and (3.14) shows that, with a 
single mass measurement, the normalization and slope of the single- field model, /nl(&p) 
and n^, are completely degenerate, and only /nl(^) at a single k value is measured. By 
adding a wide range of masses, this degeneracy is broken. This is shown in the left panel 
of Fig. 2, where narrow mass bins only constrain this function near a single k value. In 
contrast, the right panel shows that for the multi-field scenario the degeneracy is broken 
even with halos of a fixed mass, as expected from Eq. (3.18). 

Similarly, Fig. 3 contains more visual information on how the degeneracy is broken 
with multiple mass measurements. The left panel shows the constraint in the fNL{k p )-nf 

6 The careful reader will notice that these six bins in mass are an oversampling of the 30 original bins 
in mass we assumed in this interval, which are a subset of the total of 40 bins in mass we assumed in 
M = [10 13 ' 5 ,10 15 ' 5 ] /i _1 M Q . 
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Figure 3: Left panel: constraints in the /nl(&p)~^/ plane in the inflaton model. Lines show 
degeneracy directions that each of six individual mass bins suffers (these mass bins correspond to 
colored curves in Fig. 2). Right panel: Constraints in the n f -rSj" plane assuming both single-field 
and multi-field models, and marginalizing over the amplitude (term /nl(^p) = £,s{k p )£, m (kp) 2 in 
Eq. (3.22)). 

plane for the DES survey, with lines showing degeneracy directions that each of the six 
individual mass bins suffers (these mass bins correspond to colored curves in Fig. 2). The 
right panel shows the constraint in the rvfi-nj"' plane, marginalized over the amplitude 
fNh(k p ) = £,s(kp)£, m (k p ) 2 , for both DES and LSST surveys. 

Clearly, even the information from large-scale structure alone offers the possibility of 
distinguishing the origin of primordial non-Gaussianity by constraining both single-field 
and multi-field model parameters simultaneously, but the most interesting level to probe 
is running of order the spectral index, n^f ~ 0(0.04). It is not completely clear if we 
can reach that level, and there are several factors that could push the predictions above in 
either direction. First, the forecasts presented here are in some sense a best-case scenario, 
given that for simplicity we did not marginalize over the standard cosmological parameters, 
and we assumed no systematic errors in recovering the power spectra of halos, only taking 
into account the statistical uncertainties. In particular, measurements of mass of clusters 
of galaxies suffer from statistical and systematic errors that are currently at least at the 
10% level per cluster. On the other hand, constraints presented might be reached in the 
near future because we do not expect that serious degeneracies exist between the non- 
Gaussian and other cosmological parameters [15]. One exception might be the Gaussian 
bias 6g, which will need to be measured concurrently rather than predicted by theory as 
we assumed here. 

In addition, we caution that our simulation results do not agree quantitatively with 
the analytic prediction (see the next section). In fact, the simulations find a substantially 
stronger dependence on mass that what is predicted. If the simulations are proven correct, 
then the effects shown above will be easier to distinguish, which is very encouraging for 
distinguishing between different mechanisms that may lead to large local non-Gaussianity. 

Finally, the constraints above are from bias alone, but the Planck satellite will consid- 
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erably improve constraints on /tvl from the CMB through measurements of the bispectrum 
itself, so that any scale dependence may also be constrained [79]. Whereas current con- 
straints come from k ~ 10~ 3 7i Mpc, Planck constraints will extend to higher t and should 
overlap with constraints from the bias. Extant CMB analyses [79] have used a different 
parametrization of possible scale-dependence, and so it would be interesting to repeat this 
with our ansatz. As pointed out by [70] and explored in detail by [96], even more infor- 
mation can be extracted from joint constraints on models where two fields contribute to 
the curvature fluctuation. Multiple observations can separate the inherent size of the non- 
Gaussian interaction in one field (/^) from the rescaling by the fraction of power from 
the non-Gaussian field, which combine as shown in Eq. (2.16) to give the amplitude 
of the non-Gaussian term in the bias, = f^L^ 2 {k), that we have constrained here. 

4. Simulation results 

To check the dependence of the effective /nl on the tracer mass, we generated initial condi- 
tions with a non-zero bispectrum of the form shown in the first line of Eq.(3.3) (the single 
field model) with scale-dependent amplitude as defined in Eq. (3.2). This is equivalent to 
Eq. (3.4), so our function £ s (k) corresponds to a commonly used definition for /jvx(&) (see 
also Eq. (3.22)). This is the simplest possible model, but serves to check the predictions for 
how the effective coefficient of the non-Gaussian bias ("/nl") var i es with the mass of the 
halo. We have also performed a small number of simulations using the bispectrum form 
in the second line of Eq.(3.3), to verify that the bias has the expected scale-dependence. 
These simulations confirm that in such models Ab no longer simply scales as k~ 2 on large 
scales, but has n ™ dependence as well. However, at large enough values of the running 
to be distinguished by our simulations, second order effects are significant. For now we 
focus on the first scenario which is simpler and already uncovers a disagreement between 
the analytic predictions and the numerical results. 

To perform these non-Gaussian simulations, we first generated a realization of a Gaus- 
sian random field $>(x) with amplitude chosen so that = 0.8 and with spectral in- 
dex n s = 0.96. Then we squared the field, Fourier transformed and multiplied by the 
scale-dependent /nl shown in the second line of Eq. (3.22). Finally, we transformed this 
component back to real space and added it to the Gaussian piece. We evolved the result- 
ing non-Gaussian field forward from the scale factor of a = 0.005 using a flat cosmology 
consistent to WMAP7 best fit values (fi m = 0.27, h = 0.7). 

We ran 8 realizations each of Gaussian and non-Gaussian initial conditions, includ- 
ing ^ s (kp) = f NL (k P ) = 100 with n ( f s) = 0,0.6, /nl(M = 300 with n { f s) = 0, ±0.6, and 
/nl(^p) — 630 with rij = 0, —0.6. All cases had pivot point k v = 0.04 Mpc^ 1 . The box size 
was 2400 h~ l Mpc, with (1024) 3 particles, giving a mass per particle of 9.65 x W 11 h~ 1 M Q . 
Although our low-mass halos do not have many particles in these simulations, we used a 
few Lbox = 520fa _1 Mpc simulations (with fNL(k p ) = 300) where these halos were well- 
resolved to verify our results. These simulations were performed on the SciNet machines, 
where each run took about 3.5 hours on 16 nodes. 
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Figure 4: Dependence of scale-dependent non-Gaussian bias on mass, inferred from simulations. 
Left panel: Simulation results for the non-Gaussian contribution to the bias of halos with mass 
4 — 8 x 10 13 h^ 1 Mq. The black circles points have constant £ s (fc p ) = /ivx(fcp) = 300, the blue 
squares have the same £ s (fc p ) but rSf' = —0.6, and the red triangles have — 0.6. Error bars are 
sample variance across several simulations with the same parameters. Right panel: The same set of 
curves for halos with mass 32 — 64 x 10 13 /i _1 M . The scatter here is larger than in the previous 
plot since there are fewer objects at this mass. 



We find that the simulations with constant /nl are offset from the analytic expectation 
at small A; by a factor that is nearly constant with mass and is less than one for both positive 
and negative /nl- This is consistent with findings by other simulations, and the behavior 
of the offset was studied in detail by Pillepich et al. [97] and Giannantonio and Porciani 
[91]. However, as discussed above, we will effectively fit this offset out and examine only 
the difference in behavior between our /jvl constant simulations and those with running. 

From the simulations with scale-dependent non-Gaussianity, we find that different mass 
objects are indeed sensitive to an effective /nl that depends on the scale of the object and 
which increases (decreases) for positive (negative) running as the mass and size of the 
object decreases. Figure 4 illustrates this effect: the non-Gaussian term in the bias for 
small mass objects has a smaller (larger) amplitude for positive (negative) running than 
for constant /jvz, (left hand panel). The curves converge for larger mass objects (right hand 
panel). The bias correction Ab is calculated from the difference between the matter-halo 
cross correlation in a Gaussian simulation and the non-Gaussian case built from the same 
Gaussian realization, then averaged over realizations. 

The qualitative effect we expected is present, but for some halos the magnitude of the 
effect is not well predicted by the analytic expressions from Section 3. Figure 5 shows the 
deviation between simulation and prediction for /nl = 300, = ±0.6 (the same generic 
trend was seen in all parameter sets). We plot the ratio of the non-Gaussian correction 
with running to the non-Gaussian correction for constant Jnl'- 

^ sim = fe(/jVL = 300,4 8) =0.6)-b(/ jVL = 0) = Abjnf) 

b(f NL = 300, nf = 0) - b(f NL = 0) Ab(n ( f s) = 0) ' 
This is compared with the theoretical expectation calculated from the small k limit of 
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Figure 5: Simulation results for the scale-dependent non-Gaussian bias compared to theory. In the 
left panel, the vertical axis shows the mass-dependent ratio of the bias for non-Gaussianity that runs 
compared to the /jvl constant case, measured from /jvx = 300 simulations at z = 0. The upper lines 
have n>p = 0.6 and lower lines show riff' = —0.6. Redshifts z = (blue, higher values of u(M), and 
z = 1 (red, lower values of a(M) are shown. The dashed lines are the analytical prediction, showing 
that agreement is better at small a(M). The right panel shows the same information, but plotted 
as a function of mass. Now the theoretical prediction (solid black lines) is redshift independent. 



Eq.(3.14). The curves are plotted as a function of a(M). As the figure demonstrates, the 
simulation results agree well with our analytic model in the high-mass limit a(M) <C S c , 
but towards lower masses (e.g. a(M) > 0.8) the simulations produce a stronger effect than 
Eq.(3.14) would predict. Note that the discrepancy does not appear at a fixed mass, but 
rather at a fixed a(M). This is illustrated in the right panel of Figure 5, which is identical 
to the left panel of Figure except that now mass M is the abscissa. The figure shows 
that for fixed mass M, the simulations agree with Eq.(3.14) at high redshift, but begin to 
disagree at low redshift as a(M) grows. 

One very plausible explanation for this discrepancy at low mass (a > 0.8) is that the 
profiles of the peaks that produce halos begin to change as a increases. As we have argued, 
the non-Gaussian bias of halos of mass M is sensitive to the value of /jvz, at some effective 
k oc Af -1 / 3 . Implicit in this scaling is the assumption that the profiles of peaks that collapse 
into halos are similar at different masses, just rescaled in size. However, we know that this 
assumption is incorrect. Bardeen et al. [98] argued from Gaussian statistics that as a(M) 
increases, the peaks that collapse into halos generally become steeper. N-body simulations 
confirm the presence of this effect, but show that it is much stronger in magnitude than 
predicted by Bardeen et al., apparently due to environmental effects during halo formation 
[99]. Because peaks at high a(M) are much steeper than rare peaks at low a(M), they are 
sensitive to non-Gaussianity at higher wavenumbers, even at the same peak size R. For 
scale-independent /jvx, this change in peak profile has no effect, but for nonzero nf it can 
dramatically enhance the mass dependence of non-Gaussianity, as our simulations show. 
It remains to be seen whether the magnitude of the change in peak profile can account for 
the discrepancy between our simulations and Eq.(3.14); this is work in progress. 
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5. Conclusions 



In summary, we have introduced a generalization of the local ansatz, Eq. (2.9), that is 
a symmetric, factorizable function of the momenta and includes scale-dependent non- 
Gaussianity. This more general expression is motivated by natural features of models 
that give an observably large amplitude for local type non-Gaussianity, and distinguishes 
between non-Gaussian curvature fluctuations generated by a single field and multiple fields. 
If only one field contributes to the curvature fluctuation (and is different from the inflaton 
so that the non-Gaussianity may be large), the scale-dependence of the non-Gaussianity 
characterizes the self-interactions of the field. If two fields contribute to the curvature 
fluctuations, scale-dependence indicates how the ratio of power in the fields changes, which 
is a function of how different the potentials are. If local non-Gaussianity is large enough 
to be observed, such scale-dependence is as natural as running of the power spectrum. 

Models with scale-dependent local non-Gaussianity can generate two signatures in 
the non-Gaussian contribution to the halo bias. First, the non-Gaussian term may be 
proportional to an effective /nl related to the amplitude of the bispectrum on the scale of 
the object so that different mass objects have a different amplitude correction. Second, the 
k~ 2 behavior of the non-Gaussian bias can be modified to k~( 2 ~ n f ^ (where < 1), 

and one should expect the first effect to accompany this one. 

We have used N-body simulations to verify that different mass objects do indeed have a 
non-Gaussian bias proportional to an effective /nl that varies with the mass of the object. 

(s) (m) 

It is interesting that the simulations show that scale dependance with Tij < can erase 
the scale dependent effect on the bias for some range of masses, highlighting the need for 
analysis using multiple tracers of different mass. However, the quantitative result for halos 
at large cr(M) is not well predicted by our analytic expressions. We have speculated that 
the origin of this discrepancy may be related to differences in the initial peak profiles of 
the halos, but leave a detailed investigation for a later work. 

Future surveys are sure to bring interesting results. Using the analytic predictions, we 
find that they may be able to distinguish the different pieces of our generalized local ansatz, 
and so different origins of local non-Gaussianity, especially if the running is somewhat large 
|- n ( s ).( m ) ^ 0(o.l)). However, the existing analytic expressions predict a weaker effect than 
we see in the simulations, and our forecasts only account for constraints from massive 
groups and clusters of galaxies, neglecting the (potentially) greater sensitivity to running 
possible when galaxy correlations are included as well. Our forecasts for future surveys 
should therefore be taken as a lower limit on the potential to observationally distinguish 
these features. 
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